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Abstract 

We construct a class of exact ground states of three-dimensional periodic Anderson models 
(PAMs) — including the conventional PAM — on regular Bravais lattices at and above 3/4 filling, 
and discuss their physical properties. In general, the / electrons can have a (weak) dispersion, and 
the hopping and the non-local hybridization of the d and / electrons extend over the unit cell. The 
construction is performed in two steps. First the Hamiltonian is cast into positive semi-definite 
form using composite operators in combination with coupled non-linear matching conditions. This 
may be achieved in several ways, thus leading to solutions in different regions of the phase diagram. 
In a second step, a non-local product wave function in position space is constructed which allows 
one to identify various stability regions corresponding to insulating and conducting states. The 
compressibility of the insulating state is shown to diverge at the boundary of its stability regime. 
The metallic phase is a non-Fermi liquid with one dispersing and one flat band. This state is also 
an exact ground state of the conventional PAM and has the following properties: (i) it is non- 
magnetic with spin-spin correlations disappearing in the thermodynamic limit, (ii) density-density 
correlations are short-ranged, and (iii) the momentum distributions of the interacting electrons 
are analytic functions, i.e., have no discontinuities even in their derivatives. The stability regions 
of the ground states extend through a large region of parameter space, e.g., from weak to strong 
on-site interaction U. Exact itinerant, ferromagnetic ground states are found at and below 1/4 
filling. 

PACS numbers: PACS No. 71.10.Hf, 05.30.Fk, 67.40.Db, 71.10.Pm, 71.55.Jv 
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I. INTRODUCTION 



The periodic Anderson model (PAM) is the basic microscopic model for the investiga- 
tion of compounds with heavy-fermion or intermediate valence properties such as cerium 
or uranium 1 . The model describes / electrons which interact via a strong on-site Coulomb 
repulsion U and hybridize with noninteracting d electrons. In the simplest case the / elec- 
trons are assumed dispersionless, the hybridization purely local, and the d electron hopping 
nonzero only between nearest neighbor sites. However, for real systems this is an oversim- 
plification since there is experimental evidence for (i) a weak, but finite dispersion of the / 
electrons, especially in uranium compound a 2 ' 3 ' 4 , (ii) nonlocal contributions to the hybridiza- 
tion, and (iii) hopping of the d electrons beyond nearest neighbors^. 

Recently the PAM was employed to study the dramatic volume collapse at the a — > 7 
transition in cerium compounds.—^ These investigations called attention to the possibility 
of a Mott metal-insulator transition in the PAM. In fact, a remarkable similarity between 
the Hubbard model with nearest-neighbor hopping 9 - and the PAM with nearest neighbor 
hybridization and d hoppin g 9 ' 10 ' 11 ' 12 was found. These results show that the range of the 
hopping and hybridization in the PAM are quite important but still poorly understood. 

In this situation exact results on the existence of insulating and metallic phases in the 
PAM and their dependence on the hopping, hybridization, and interaction parameters are 
particularly desirable - especially in three dimensions. So far exact results for the PAM 
were mostly limited to special regions of parameter space, namely, for infinite repulsion of 
the / electron a 13 ' 14 ' 15 , and for finite repulsion in low dimensions D = 1, 2 16 i 17 i 18 i 19 

In this paper we not only present details of the construction and the physical properties 
of a class of exact ground states of three-dimensional (D = 3) periodic Anderson models 
reported in Ref.— , but extend the range of applicability of our approach substantially. In 
particular, we (i) demonstrate the uniqueness of the metallic and insulating solutions discov- 
ered at 3/4 filling 2 ^, (ii) explicitly present and analyze the non-linear matching conditions 
connecting the starting Hamiltonian to the transformed Hamiltonian, (iii) deduce the cur- 
rent operator and the sum rule for the charge conductivity, (iv) derive several local and 
global expectation values such as the magnetization of the system, (v) calculate correlation 
functions, (vi) extend the solutions to the conventional PAM case, and (vii) show that, by 
employing different procedures to cast the Hamiltonian into positive semi-definite form, one 
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arrives at exact ground states in different regions of the parameter space. 

The paper is structured as follows. In Section II we present the Hamiltonian, discuss 
its transformation into positive semi-definite form, and construct a class of exact ground 
states. Section III describes the localized solution, and Section IV characterizes the metallic 
non-Fermi liquid state. In Section V the approach is generalized, leading to solutions in 
other regions of parameter space, and in Section VI the results are summarized. Technical 
details are discussed in Appendices A - E. 



II. TRANSFORMATION OF THE HAMILTONIAN AND CONSTRUCTION OF 
EXACT GROUND STATES 

A. General form of the periodic Anderson model 

We consider a general form of the periodic Anderson model (PAM) describing non- 
interacting d electrons which hybridize with interacting / electrons. In contrast to the 
conventional PAM we do not assume the / electrons to be localized, i.e., the Hamiltonian 
is given by 

H = H + H v , (l.a) 

h = £[(4< + 4pL) + (vju^ + v:ftj^)i (i.b) 

k,<r 

Hu = Uj2h(Xl- ( Lc ) 

i 

We denote the two types of electrons by b — d, f, i.e., b\. a creates a b(= d, f) electron with 
momentum k and spin a. The corresponding particle number operators are n^ a = b^J)^, 
and the dispersion relations of the 6-electrons are given by e£. Furthermore, the hybridization 
amplitude and the local (Hubbard) interaction are denoted by Vk and U, respectively. 

In real space the model is defined on a general Bravais lattice in D = 3 dimensions, with 
a unit cell I defined by the primitive vectors {x T }, r = 1,2,3. The non-interacting part of 
the Hamiltonian, H Q , reads 

^0 = £{£[(^Vi+ r ,* + t f JlJi + r,a) + (V r dJ dlJi + r,a + flj^) + H.C.] 

+ (VodlJ^ + H.c.) + E f n( a }, (2) 
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where t b r characterizes the hopping of b electrons between sites i and i + r, V b,b ' is the 
hybridization of b and b' electrons at sites i and i + r, Vq is the on-site hybridization, and 
Ef is the local on-site /-electron energy. The separation between a site i and its neighbors 
in the unit cell is denoted by r, with r ^ 0. While the amplitudes t b are real, Vq, V b,b ' 
can, in principle, be complex (whether V b,b ' is real, imaginary or complex, depends on 
the linear combination of the corresponding electronic orbitals and hence on the lattice 
symmetry^ i 18 i 21 i 22 i 23 ) and obey the relations 

4 = E f S bJ + E(^e" ikr + ^V ikr ), V k = V + £(K d ' J f e" ikr + V r f ' d *e +ikl ). (3) 

r r 

In particular, for t{ = 0, the /-electrons are localized, and the model reduces to the 
conventional PAM. 

B. Transformation of the Hamiltonian. 

1. Representation of sites in a unit cell 

The separation from a site i in (J2J) is indicated by the vector r which corresponds to 
neighboring sites located in different coordination spheres. In our investigation r may extend 
over a unit cell I of a general Bravais lattice in D = 3. This implies 26 different inter-site 
hopping and non-local hybridization amplitudes. To avoid multiple counting of contributions 
by the H.c. term in (J2J) the vector r must be properly defined. To this end the sites within 
Ii, the unit cell defined at site i, are denoted by 17. = i + r Qj g 7 , with r a/ g 7 = ax x + /5x 2 + 7x3; 
a, [3, 7 = 0,1. As shown in Fig. 1 the eight sites r/. can be numbered by the indices 
n(a,/3,7) = 1 + a + 3(3 + Aj — 2af3 without reference to I\. Then r = r^/y — r Q/37 , with 
n(af, /3f, 7/) > n(a,/?,7), connects any two sites within a unit cell. It corresponds to half of 
the 26 possibilities, i.e., to the 13 possibilities xi,X2,X3,x 2 ±xi,X3±xi,X3±X2,X3±x 2 ±xi. 
The remaining (negative) values of r are taken into account in Q by the H.c. contributions. 

2. Transformation of H into positive semi- definite form 

To construct exact ground states the Hamiltonian H needs to be rewritten in terms of 
positive semi-definite operators. This is made possible 24 by the construction of two new 
operators - one for the transformation of ifo an d one for Hjj. 
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The first of these operators is the " unit cell operator" A\ a which represents a superpo- 
sition of fermionic operators creating d or / electrons with spin a inside every unit cell I\ 
as 

8 

Azi,<7 = [ a n,d^i+r Qj/3i7 ,o- + a nJ f\+r a ^ rj ,o\ 

n(a,(3, 7)=1 

= ( a l,dd\,a + a 2,d^i+xi,<r + a 3,d C 4+xi+x 2 ,<T + a 4,d^i+x 2 ,a + •"■ + a 8,d^i+x 2 +X3,(r) 

+ ( a i,//i[cr + a 2,//i+Xl,CT + a 3,//i+Xi+X 2 ,cr + a 4,///+X 2 ,cr + •"■ + a 8,//i+X2+X 3 ,cr) • (4) 

Because of the translational symmetry of the lattice, the numerical prefactors a* b , n = 
n(a, P, 7) are the same in every unit cell. The composite operators A\ a do not obey 
canonical anticommutation rules, since {Aj tU , A\, a ,} ^ for all I 7^ I'. This is because 
Aj a creates electrons also at the boundaries of the unit cell Ji with neighboring unit cells. 
It should be stressed that for this reason A\. a has a genuine dependence on the lattice 
structure and thereby on the spatial dimension. Furthermore the relations {A\ a , A\, a ,} = 
{A Ij(T , A ry } = and A\ <a ] = K d + K f , K b = E„=i K,b\ 2 im P l Y 

(A-U 2 = 0, (5.a) 
-i} )( A* = A I;a A\ ja - (K d + K f ). (5.b) 



The second operator, P, originates from the relation Ei = 
where 

p = EA, Pi = - (™£t + Ki) + !> (6) 

i 

and N\ is the number of lattice sites. The local operators Pi are positive semi-definite and 
assume their lowest eigenvalue (=0) whenever there is at least one / electron on site i. By 
contrast, the interaction operator Hjj itself, which is also positive semi-definite for U > 0, 
assumes the eigenvalue only if the double occupancy is exactly zero. For this reason P is 
more useful for our investigation than Hu. 

Taking into account periodic boundary conditions and allowing r to take only the values 
discussed above, (|l.al EJ) may be written as 

H = P A + UP + E g , (7) 
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where P A = A Iita A\. a , E g = K d N + UN A - 2N A (2K d — Ef), and N is the total number 
of particles which is assumed to be fixed. 

For ((7J) to reproduce (jl.aj) the prefactors a* b in A} ;tT must be expressed in terms of the 
microscopic parameters tf, t{., V r , V* , Vo, V Q * , Ef, U, for all r £ Ji, taking into account 
periodic boundary conditions. This leads to 55 coupled, non-linear matching conditions 
which can be written in compact notation as^ 



1 3 

(n. D ffim) a n+ ) b a n- ) b l 
Pi,j3 2 ,l3 3 =-1 i=l 

These matching conditions have the explicit form 



rpb,V 



j: 



6,6' 



J: 



J: 



j: 



j: 



0*l£ a 2,b' + al : b a 3,b' + O-l^G^' + a 8,6 a 7,6', ^xf+xi — °l,6 a 3,6' + a 5,6 a 7,6', 

* * * fob'** 

a l,6 a 4,6' + a 2,b a 3,6' + a 6,6 a 7,6' + a 5,6 a 8,6' , ^x 3 +xi = a l,6 a 6,6' + a^b a 7,b', 

a l,6 a 5,6' + a 2,b a 6,b' + a 3,6 a 7,6' + a 4,& a 8,&') ^x 3 6 +x 2 = a l,6 a 8,6' + a 2,6 a 7,6' ; 



Xl 



a 2 6 a 4,6' + a , 



6.6 



^xf-xx — a 2,6 a 5,6' + «3,6 a 8,6' , ^x^-xa — a 4,fe a 5,6' + «3,6 a 6,l 



* t6,6' _ * jO,o _ * 

X2+xi — a l,6 a 7,6', ^x 3 +x 2 -xi — a 2,b a 8,b', ^x 3 -x 2 +xi — a 4,6 a 6,6' 

8 

a 3,b a 5,6', Vo = - 51 CL*n,d a n,f, U + Ef = K d ~Kf, 



jb,b 
^X 3 -X 2 -Xi 



T 6,6' 



(8) 



n=l 



where J r 6 - 6 ' = -[^ + (1 - 6, = d, f. 

Details of this transformation are presented in Appendix A. 



C. Construction of exact ground states 

Apart from the constant term E g in J7J) H is a positive semi-definite operator. A state 
\ty g ) which fulfills the conditions 

A|* fl ) = 0, (9.a) 
4j*,> = (9.b) 

for all i, and which contains all linearly independent states with properties ()9.al9.b|) . will 
then be the exact ground state of H with energy E g . Since the kernel of an arbitrary operator 
O, ker(0), is defined by the linearly independent states \4>) satisfying 0\4>) = 0, the relations 
()9.a|) and (|9.b|) define the kernel of the operators P and Pa, respectively. Consequently, \^f g ) 
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spans the common part of ker(P) and her (Pa) denoted by the Hilbert space 

H g = ker(P A )nker(P). (10) 

Using this definition it is ensured that \ty g ) is the complete ground state, and that supple- 
mentary degeneracies of E g do not occur. 

Using © and O it follows that ker(P) is defined by states F*\Q) = llfJi ^ f |0), where 
Fi = (/ii,t/i[| + A*u/n) an d Aii,o- are arbitrary coefficients. Obviously F^ creates one / 
electron on every site i. Furthermore, ()5.a|) and (|9.b|) imply that /cer(P^) is defined by states 
G^|0), where & = l\i!ja(A\.^ A\. ^) creates at most two (d or /) electrons on i. Since & also 
creates contributions without / electrons, implies the (unnormalized) ground state 

\* g ) = &Fi\o) = l[[A\ iA A\ i J}}\o) (ii) 

i=l 

at = 3N^, e.g., 3/4 filling. Clearly, \^ g ) has the desired property H\^ g ) = E g \ty g ) and 
spans 7i g at 3/4 filling. Thus it is the exact ground state of H with energy E g . 

Since the operator F^ was introduced into |\l/ g ) to take into account the operator P in 
(J7J), i.e., the Hubbard interaction U, \^ g ) can only be a ground state for U > 0. 

Eq. (fTT]) implies that (i) the linearly independent basis vectors of 7i g have the form 
|^ g {a;}) = GtF^jO), where F ] {a . } = ]jf=i f!^, and (ii) G t does not contribute to the total 
spin of the ground state. The overall degeneracy and total spin S G [0, of \^ g ) are then 

determined only by the (arbitrary) set of coefficients \i\ iCJ (see also Sect. III. A. 2., Sect. IV. E., 
and footnote p£]). Consequently, the degeneracy of E g is determined by the (high) spin 
degeneracy of \^ g ). Since in \ty g ) all possible values S and orientations S occur the ground 
state is globally non-magnetic^. 

Exact ground states can also be constructed away from 3/4 filling. For example, the 
operator = YlfLi EiI A i J2b=dj,a a l,j,(M a\ with numerical coefficients ay , creates M < N\ 
additional particles into the system such that 

N A 

\%) = IM, Ai^i^iio) (12) 

i=l 

is a ground state for U > and N > 
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III. EXACT LOCALIZED GROUND STATE 



The physical properties of \^ g ) depend on the values of the coefficients a n ^ in (HJ) which 
are solutions of (jHJ) for given microscopic parameters. We now identify different solutions 
for nonlocal hybridization amplitudes with V® = Vj d = V r and discuss their physical 
properties. We start with a localized ground state at 3/4 filling. 



A. Derivation of the localized ground state 

From it follows that if a* nd and a* ^ are proportional, i.e., a* d = pa* n j for all n, the 
operators A\ a take the form 

8 

A\ u a = H a n{aAl)J E Ur at0 ,^ ( 13 ) 
n(a,/3, 7)=1 

where E\ a = (pd\ a + // CT ). The ground state \^ g ), (JHJ, then transforms into |\E% C ) = 
n^i^tAa^)! )- Evaluating the product E^E^E?, one finds 

N A 

\*ioc) = nt£/v(p4Ait + flJlAvW)- (14) 

i=l cr 

Since (^i oc \^i oc ) = (1 + \p\ 2 ) Na IISUKtI 2 + l/^ul 2 ) ^ the ground state is well-defined. 



1. The insulating nature of the ground state 

The state |\l// oc ) has exactly three particles on each site, corresponding to a uniform 
electron distribution in the system. Indeed, for hi = J2b=dj J2a n\ a one finds ni|\E% c ) = 
3\^i oc ). Denoting ground state expectation values in terms of |\E%c) by (....), one obtains 

(4A.') = 0, (d{J jy ) = 0, (ftAy) = 0, (fUiy) = 0, (15) 

for all cr, a' and all i ^ j. Hence hopping or non-local hybridization does not occur. 

By separating the Hamiltonian H into an itinerant part Hum = S r -^iim( r ) an d a com- 
plementary localized part Hi oc — H — H it i n (see Appendix B), and using (JEJ), one finds 
(Hiti n (r)) = for all r, and (Hi oc ) = E g . This clearly demonstrates the localized nature 
of the ground state. Furthermore, from (|B3|) . the sum rule for the charge conductivity 
is obtained as / °° dujRea T T {uj) = 0. Since Rea TtT {uj) is non-negative this relation implies 

8 



Rea T)T {uj) = 0. In particular, Rea TtT (0) = 0, the DC-conductivity, is also zero. The ground 
state (fTlj) is therefore insulating2£i2£ It should be noted that the nature of this state is quite 
nontrivial, since the localization of the electrons is due to a subtle quantum mechanical 
interference between states with two d and one / electron {d\ +d\ *f<[ a ) and two / and one d 
electron (/jL/J, d\ ) on every site i. 

A state with Rea TjT (uj) = for all u appears to be rather unphysical since it implies that 
not only the DC-conductivity but even the dynamic conductivity vanishes for all excitation 
energies. It should be stressed, however, that the relation Rea TiT (uj) = was derived in 
the framework of the Kubo formula for the charge conductivity, i.e., within linear response 
theory. Consequently, this result is not valid at high excitation energies uj. 



2. Global magnetic properties 

The expectation value of the spin^ in terms of the ground state (j!4j) in Cartesian coor- 
dinates is found as 

{) ti^M 2 + M) k 2(K T p + | m H + h*(M + MY K 

where |x| = |y| = |z| = 1. The total spin is seen to depend on the arbitrary coefficients 
Here the site dependence of the ^i- ha coefficients should be stressed. Namely, by 
choosing /x i (T = fi a one obtains 

~ 2 iv A jv a /A g v _iy A (N 2 -N a ) n7 x 
(S) "^ ( ^ + ) ' <5> ~1"(N 2 + N 2 )' { ) 

which represents a ferromagnetic state with maximal total spin S/N^ = 1/2, leading to 



(S 2 )/Na = 1/2 in the thermodynamic limit; is seen to influence only the orientation of 
S. 

The minimal total spin can be found by considering two distinct subsystems of arbitrary 
shape, both containing the same number of lattice sites Taking /i i (T = n a ,iii_ a = 

in each subsystem one finds 



lS " ; ' (S x ) = (S y ) = (S z ) = 0, (18) 



which implies zero total spin (i.e., a global singlet state) in the thermodynamic limit. De- 
pending on the choice of the parameters {/Ui jCT } all values of S between these two extreme 
values for S, and all orientations of S, can be constructed (see Sec. II. C). 



3. Local magnetic properties 



Analyzing the local magnetic properties of the ground state one observes that the expec- 
tation value of the double occupancy per site for both b = d, f electrons is smaller than unity 
since {n\+n\\) = (5bj + \p\ 2 Sb t d)/ (1 + \p\ 2 )- Consequently, each site carries a local moment. 
Indeed, irrespective of the values of fi- ha and p one has (S 2 ) = 3/4 on each site i; this is the 
result of a quantum mechanical superposition of the corresponding contributions of d and / 
electrons. The / and d moments do not compensate each other locally since 

((s^) + (s M »(£ |/v| 2 ) = x^uKt) + yM/*uKt) + z lRt|2 o — » ( 19 ) 

where (S^d) = |p| 2 (Si,/) holds. Furthermore, taking into account fixed (but arbitrary) fi- ha 
the spin-spin correlation function for i ^ j is found as 

An average over all possible values of /Xi )CT , /ij jCr therefore implies (SiSj) = 0. Therefore, 
in spite of the existence of local moments the system is globally non-magnetic; this is a 
consequence of the large spin degeneracy of the ground state. 
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B. Solutions of the matching conditions 

The matching conditions (JBJ) for the nonlocal hybridization amplitudes along the space 
diagonals of the unit cell read 

— ^x 3 +x 2 +xi = a *,d a 7J = a l,f a 7,d: ~ ^x 3 +X2-xi = a 2,d a S,f = a 2,f a S,d, 

— V r x 3 -x 2 +xi = fl4,d a 6,/ — a 4,/ a 6,d, — Kc3-X2-xi = O^OS,/ = a 3J a 5,d- (21) 

They hold for all {x T }, r = 1,2,3 and hence imply a* n Ja* n j = a n ' j( i/a n >j = p = p*. Since 
for real p, (JBJ) leads to the relations 



*/ _ £ 1/ - £ y + _ iz£ 

p z p Vq p 



(22) 



real p are seen to imply real hybridization amplitudes. In addition to (|22|). (JSJ) yields the 
following system of 14 coupled nonlinear equations 

— * X1 = a l,d a 2,d + a 8,d a 7,d + 0>Xd a 3,d + a l,d a 6,d' 
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— ^x 2 ~~ a *l,d a ^4 + a 6 t d a 7,d + a 2,d a 3,d + a 5,d a 8,d, 

— £ Xg = a* d a5 !( i + a* zd a 7 ^ + a^ d a 6i d + a^ag.d, 

— r x 2 +xi = a l,d a 3,d + a 5,d a 7,d, — r x 2 -xi = a 2,d a i,d + a 6,d a 8,c(, 

~^x 3 +xi = a l,d a 6,d + a l,d a 7,d' _i x 3 -xi = a 2,d a 5,d + 0* 3d a 8 ,d, 

x 3 +x 2 = a l,cZ a 8,d + a 2,d a 7,rf; — r x 3 -x 2 = a 4,d a 5,d + a 3,d a 6,d, 

— ^x 3 +x 2 +xi = a l,d a T,di ~ ^x 3 +x 2 -xi = a 2,d a 8,d, ~ ^x 3 -x 2 +xi = a 4,d a 6,di 

8 

-*x 3 -x 2 - Xl = a s,d a 5,d, pVb = - 2 l a ",d| 2 - ( 23 ) 

n=l 

They determine the unknown complex coefficients (i.e., 16 unknown real values) from 
the input parameters tf. A study of the possible solutions shows that for \p\ > 1 the 
relative sizes of the hopping and hybridization amplitudes are physically very reasonable, 
e.g., < |t£ |, |txi+x 2 l < l^xi I; l^xi+x 2 l < l*xj- That is, they decrease with distance, 
and the magnitude of the amplitudes of the d electrons is larger than those of the almost 
localized / electrons. 

Based on (J7J), the corresponding ground state energy becomes E g /N^ = —U + (1 — 
2/p 2 ) J2n=i \a n ,d\ 2 - Depending on the solution, E g has a nontrivial structure which will be 
analyzed below. 



1. Solution for a cubic lattice 



For a simple cubic lattice one has t Xl = t X2 = t^ 3 = tf, t X2±Xl = t x , s±yii = t X3±X2 = t%, 



£x3±x2±x! = Then has a solution 



a i,d — di, a 2,d — o^d — a 5 ^ d — ua 1} a 3d — — a s,d — u a x , a 7d — u ai, (24) 

where a\ is an arbitrary, real quantity which is fixed by the energy unit. Furthermore, u 
is real and is determined by the parameters entering in H through (J23I24)) . Due to the 
almost localized nature of the / electrons their nearest-neighbor hopping amplitude can be 
expected to be much smaller than that of the d electrons. Indeed, for z = |t{/txl < 1/2 the 
ground state energy is given by 
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where u = (1 - y/1 - Ay' 2 )/(2y), and y = \t%/t(\ G (0, 1/2] holds. We see that |^ Zoc ) ceases 
to be the ground state for y = |tf/tf| > 1/2 = y c . This corresponds to rather strong 
next-nearest neighbor hopping of d electrons. Apparently, at y c = 1/2 a different — most 
probably itinerant — phase becomes stable. At y c the ground state energy E g (y) has a 
finite value, but its derivative diverges due to du/dy = +oo (see Fig. 2). Since the size of the 
hopping element may be tuned by pressure, the infinite slope of E g at y = y c is expected 
to correspond to an anomaly in the volume, or the compressibility, at a critical pressure P c . 
Such a feature is indeed observed in some heavy-fermion materials^. 

2. Solutions for non-cubic lattices 

Similar results may be deduced for other lattice structures. In the most general case, 
i.e., when all hopping amplitudes are different, the ground state energy for the localized 
solution becomes E g /N\ = —U + (1 — 2|t£ \/\t^ |) Zln=i l a n,d| 2 - When at least one of the 
terms |a n) d| 2 (see flC2|) in Appendix C) is mathematically no longer defined the localized 
solution becomes unstable. Except for accidental cancellations in the ground state energy 
an infinite slope of E g as a function of the hopping amplitudes is also found in this more 
general case (see Appendix C). 

IV. EXACT ITINERANT GROUND STATES 

The localized ground state discussed above has exactly three electrons per site. In general, 
however, the inter-site hopping and hybridization will lead to a variable number of electrons 
on each site. In that case the ground state \^f g ), (fTTj) . becomes conducting. We will now 
describe solutions of this kind. 

A. Solution for a n> & ^ for all n, b. 

To solve the matching conditions (JBJ) for the case where a n ^ ^ for all n, b, we define 

a* 

Pn = -r~,n = 1,2,3,4, (26. a) 

a nJ 

p n , = ^,n' = 5,6,7,8, (26.b) 

a n',f 
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and consider again = V/ ,d = V r . Eq. (|26.b|) is seen to contain coefficients a n fi instead 
of a* b since (|2*Tj) holds as well. An itinerant solution is obtained by choosing p n = p = —p*, 
i.e., imaginary p. For this choice (jSJ) leads to 

£ = fjli. ^ 2 = -^, K 2 =<, ^=0, (27) 
where I - ! = x 1; x 2 , x 2 ± x 1; r 2 = x 3 , x 3 ± x 2 , x 3 ± x 1; x 3 ± x 2 ± x 1; and V T1 = follows from 
V ri = —V ri . As discussed earlier, imaginary p imply imaginary V r . The local parameters U 
and Vq become 

U 2 -l 8 8 4 

[J + E^^^^la^l 2 , = E l«Ml 2 - E l«Ml 2 - (28) 

\P\ n=l n=5 n=l 

The remaining relations following from (|SJ) are 

— t Xl = a* ld (X24 + a 8,d a 7,d + «4 j( ifl3,d + a 5,d a 6,d; 
— ^x 2 = a l,d a 4,d + a 6,d a 7,d + 0-2,d a 3,d + «5,d a 8,d! 



-£ Xg — a* d a5 j( i + a* zd <XT y d + a 2 da 6iC 2 + 04,. 



— r x 2 +xi = a l,d a 3,c( + a <o,d a 7,di ~ r x 2 -xi = a 2,d a i,d + a 6, 
— ^x 3 +xi = a l,d a 6,d + «4,d a 7,d ; ~4 3 -xi = a 2,d a 5,d + «3,d a 8,d; 
— ^x 3 +x 2 = a l,ci a 8,rf + a 2,d a 7,di ~ ^x 3 -x 2 = a l,d a 5,d + 0*Z4 a &,di 
~ ^x 3 +X2+xi = a l,d a 7,di ~ ^x 3 +x 2 -xi = a 2,d a 8,d; 

— ^x 3 -x 2 +xi = a 4,d a 6,d, — ^x 3 -x 2 -xi = a 3,d a 5,d- (29) 

Furthermore, = implies 

Ol )C i a 2,(i + a 4,<i a 3,d = a 5,d a 6,d + a 8,d a 7,ci, a l,d a 3,<i = a 5,d a 7,d) 

a l,d a i,d + a 2,d a 3,d = a 6,d a 7,d + al,d a 8,di a *2,d a ^d = «6,(i a 8,d- (30) 

From (jH IHl) it follows that 

a i,d = °1j a 2,d = a 3,d = u 2 a 1 , a Ad = ua 1 , 

a* 5d = ua\, a* &4 = u 2 a\, a 7:d = ua x , a Sid = oi, (31) 

where w is real, |w| 7^ 1, Vo = 0, and X] 8 =1 |a n ,d| 2 = 2|ax| 2 (l + u 2 ) 2 . These relations can 
only be fulfilled for anisotropic hopping and hybridization amplitudes, e.g., for vanishing 



hybridization in the basal (x, y) plane. Namely, (|29H31|) yield \u\ = y |£x3-x 2 Ax3+x 2 l an d 
Kx T l = l*x / I- The anisotropy in the hopping amplitudes is seen to start at the level of 
next-nearest neighbors. The stability region of this phase is presented in Fig. 3. 
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B. Exact itinerant ground state of the conventional PAM. 



The solutions obtained sofar, namely (|21H25p . and (|27H31|) . require the / electrons to be 
itinerant, i.e., t( ^ 0. We will now show that exact itinerant ground states can even be con- 
structed for the conventional PAM, i.e., for t{ = 0. This requires nonzero <i-electron hopping 
up to next-nearest neighbors together with local and nearest-neighbor hybridizations. 

Such a solution can be constructed, for example, for nonzero coefficients 
a i,f, a i,d, a>2,d, G4,<2, a>5 t d, with the remaining coefficients a n) & = 0. This is realized in a tilted 
unit cell (see Fig. 4) where the distances between lattice sites with indices (n = l,n = 

3) , (n = 1, n = 6), (n = 1, n = 8) and corresponding hopping elements (t X2+Xl , £ X3+Xl , ^x 3 +x 2 ) 
are considerably larger than the distances between lattice sites with indices (n = 2, n — 

4) , (n — 2, n — 5), (n — 4, n — 5) and hopping elements (t X2 _ Xl , £ X3 _ Xl , *x 3 -x 2 )- For this 
reason the former hopping elements are neglected (i.e., put to zero). In this case (JHj) reduces 
to the following eleven equations 

C Xl = ~ a l > d a 2,di r x2 = ~ a l,d a <l,d' ^x 3 = ~~ a l,d a 5,dj 

*x 2 -xi = ~ a 2,d a ^di *x 3 -xi = ~~ a 2,d a 5,di *x 3 -x 2 = ~ a 4,d a M> 

= -</«2,d, V^* = -a^a^, V^/ = -a* )/ a 5)(i , 
\Z = -al d a 1>f , U + E f = K d -\a ljf \ 2 , (32) 

where = \ai id \ 2 + |a 2> rf| 2 + |a 4irf | 2 + |a 5j d| 2 , and all other amplitudes are identically zero. 

Considering only the simplest case, i.e., tf = t^ T , V\ = Vjf, r = 1, 2, 3, and t\ = t Xr _ x ,, 
r > r r , with real hopping amplitudes^, the corresponding stability region corresponds to 
the surface in parameter space (see Fig. 5) described by 

^T^z + VA (33) 

\ti\ X 

where x = l^l/l^il an d D — l^il/l^il- Further properties of the solution are discussed in 
footnotes [ 31 i 32 ]. The fact that a* n d = pa* n j with p = p* does not hold for all n = 1,2,.., 8, 
implies a variable number of electrons at each site, i.e., an itinerant ground state. 



C. Diagonalization of the Hamiltonian 

To describe properties of the itinerant state, a k-representation is more suitable. This 
can be introduced without restriction on the a n f, coefficients. Therefore the two conducting 
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solutions presented in Sections IV. A. and IV. B. can be treated simultaneously. 
Denoting the Fourier transforms of A\. a by A} Va = J2b=dj a k6^L-> with 

a k,6 = a l,b + a 2,b e 1 + a 3,b e 1+ 2 ^ + d^jfi 2 + dr )b e 3 

+ a* jfe e ik(xi+X3) + a; ib e ik(xi+X2+X3) + a* jfe e ik(x2+X3) , (34) 
one may define new canonical Fermi operators describing composite fermions Cs^a, 5 = 1,2, 

{C^ker, CV.k'o-'} = ^<5,5'^k,k'^o-,cr', {C^ka, CV,k'<x' } = 0. Hero 32 ' 33 

C'l.ker = \/^k^4ko- = \/^k( a k,d^kcr + Ctkj/kcr), C^kcr = \[Rk{o J \ L jd] ia — a^dfva) ■, (35) 

where i?^ 1 = X^=d,/ |ok,b| 2 - Then the Hamiltonian becomes^ 

H = H 9 + UP- UN A , H g = Y\{K d - R^C^Am + K d Cl M C 2 ,U (36) 

k,cr 

In the ground state where P\^f g ) = the Hamiltonian H reduces to H g . Hence the composite 
fermion operators introduced above indeed diagonalize H . There are two bands, the lower 
one having a dispersion £ lk = K d — R^ 1 , while the upper one is dispersionless, i.e., flat, 
since ^2,k = K d =const. The lower (upper) diagonalized band contains fermions created by 
(C^kff)) respectively. Since the total band filling is 3/4 the lower band is completely 
filled, while the upper band is half filled. Thus the Fermi energy is given by Ep = Kd- 
Band structures containing a partially filled flat band around Ep have indeed been observed 
experimentall y 3 ^ 3 ^ 3 ? . 



D. Conductivity of the itinerant state 

The ground state expectation value of Hum, the itinerant part of H, is obtained as 

(Hu m ) = - E Kd \ a ^ + K ^/ < 0. (37) 



This is in contrast to the localized case where (Hiu n ) = 0.— 

Since (H itin (r)) ^ the charge sum rule (jB3|) implies J °° duRea TiT (u) ^ 0, i.e., the 
(dynamic) conductivity is in general non-zero. However, the sum rule does not allow us to 
draw conclusions about the DC- conductivity cr(0). This becomes possible if we calculate the 
chemical potential of the system as a function of the particle number. To this end we observe 
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that in the case of a variable number of particles per site the itinerant ground state is defined 
at and above 3/4 filling, such that \^f g ) in (jllj) can be generalized to fillings beyond 3/4 via 
(112)) . Using (J7j) this allows one to calculate the energy E g for different particle numbers. In 
particular, one finds fi + = E g (N + 2) - E g (N + 1) = K d , fi~ = E g {N + 1) - E g (N) = K d , 
i.e., yU + — fi~ = 0. Therefore, the described itinerant solutions are indeed conducting^. 



E. Magnetic properties. 

Using k-space notation the unnormalized ground state, (jllj) . can be written as^S 

N A N A N A 

« = [n4 T 4j{IlE eiki (^T/;T +A*u/li)]}|0>, (38) 

k i k 

where the sum and product over k extend over the first Brillouin zone, and the set {/ii >(T } is 
arbitrary. Using ([3*8)1 we will now calculate ground state expectation values of the spin for 
different sets of {ii- ha }^ 



1. Maximum total spin 

As in the localized case = /v corresponds to maximum total spin. Defining F^ b = 
(A*T^kT + AH^kj.) an d employing (F^. b ) 2 = the product over sites i in (j"~J"J) may be written as 
nf A [Sk^ e * kl -^k,/] = ^n^-^k/' where Z is defined in footnote (^). Then the normalized 
ground state becomes 

9 kV(iMTl 2 + i^| 2 )(l«M| 2 + |a kj /| 2 ) 
The spin expectation value then follows as 

2(S) = mtf + mni ij^fii - Mj^j M 2 - N 2 Uq) 
Na N 2 + N 2 y N 2 + M 2 kl 2 + N 2 )' 1 J 

resulting in |(S)/iY A | 2 = 1/4, and (S 2 ) = (iV A /2)(iV A /2+l). Thus we recover (JTIJ), the results 



for the localized solution. The maximum total spin is again given by y (S 2 )/iV A = 1/2 in 
the thermodynamic limfb 2 ^. 
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2. Minimum total spin 



To determine the minimum spin value we proceed as in the localized case, i.e., divide the 
system in two sublattices X>j and T>^ containing both lattice sites, such that for all 

i n G T>i and j n G T>[ we have j„ = i n + R, where R is a fixed Bravais vector- 41 . Choosing 
A*i,CT = A*i,-(7 = for in ([5S]) the unnormalized ground state becomes 

A/a jv a 

N A -f- N A 

\%) = ni4An n (E^ k v kT )][ n (E^%)]io>- (4i) 

k ieDi k jeD 2 k 

Details of the calculation of ground state expectation values in terms of (14 lj) are presented 
in Appendix D. Based on flDHJ) one finds (S) = 0, and using (fDljl . yftp)/N A < v / 3/(2 v / A 7 a). 
As a consequence, this itinerant ground state has zero total spin in the thermodynamic limit. 

Between the two limiting cases of the total spin discussed above all values of S may exist, 
depending on the choice of the set The statement regarding the degeneracy of the 

ground state presented below (JTBJ) holds in the itinerant case as well. 

F. Momentum distribution function 

To calculate ground state expectation values of k-dependent operators involving Cs,ko- 
operators, (jSHJ), we need to express \^ g ), (jSHJ), in terms of Cg ka , too. For the unnormalized 
ground state one finds 

N A N A N A 

w = [nc , ik T c , i t ,k i ]{n[E^^ ki ( / . i , T c' 2 t ikT + /i U ci ki )]}io>, (42) 

k i k 

where {yUi, CT } is arbitrary, and X k = — a^ aVRk', in the following we consider X k ^ for all 
ki2. From (jUJ) it follows that 

dl^Cx^lVg) = 5^5^ g ) , Cl ka C 2 ,v*>\y 9 ) = 0- (43) 

Details of the calculation of ground state expectation values of C\ ka C2,v.a operators are 
presented in Appendix E. Using (|4H|E8jl . one finds in particular 

(CI m Csm) = Ks + \hs (44) 

for all k. This expresses the fact that the lower band (5 = 1) is completely filled, while the 
upper band (5 = 2) is flat and half filled. 
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Employing flH| the momentum distribution of b = d, f electrons is then obtained as 3 ^ 

<« = (AX«) = \ + \R*WA\ (45) 

which implies n ka = J2b=dj n ka = 3/2. Since the coefficients -Rk|ak,b| 2 = |ak,fe| 2 /(l a k,d| 2 + 
|ak,/| 2 ) are regular functions of k, this is also the case for n k(T and (see Fig.6) 43 . Conse- 
quently, the momentum distributions of the electrons in the interacting ground state have no 
discontinuities. Since the ground state is nonmagnetic and metallic, the system is therefore 
a non-Fermi liquid. This is a consequence of the macroscopic degeneracy of the electrons in 
the upper band. Due to the flatness of the upper half-filled band all k states are equivalent. 
Hence, even if a Fermi energy Ep = exists the Fermi surface, and the Fermi momen- 
tum are not defined. Consequently, Luttinger's theorem^ 4 - does not apply. Non-Fermi liquid 
properties connected to a flat band have also been observed in other investigation o 4 ^ 4 ^ 4 ? . 

G. Correlation functions 

To characterize the itinerant ground states further we now calculate correlation functions 
of the interacting systems in D = 3. This is made possible by the explicit form of the exact 
ground states. 

1. Density- density correlation function 

Using (|4^j) the density-density correlation function 

Pn,n(r) = ^^((nin i+r ) - (hi)(n i+r )), (46) 

where h\ = J2o-12b=d,f mSi J be written as 

Pn,n(r) = e l(kl ~ k2)r yj R kl i? ka -Rk 3 -Rk 4 [(^kg ,d a k 4 ,d + a k 3 ,/ a k 4 ,/) X 

A a, a' ki,k2,k3 

( a ki,d a k 2 ,rf + a kiJ a k2j)(( ( ^2,k 3( 7'C'2,k4 ( T'C'2,kiaC'2,k2 ( T) _ (^ ! 2,k 3 a'^2,k4,^) (^.kja^Ucaff)) + 

(«k 3 ,d a k 4 ,/ - a k 3 ,/a k4i d)(a ki / a k2)d - a^ i4 a^ 2 j)(Cl Maa ,C2^ ia 'Cl Ma C ltk2a .)}, (47) 

where k.4 = k3 + ki — k2. Employing (|E6|E12J) in the first term and (j43JEtij) in the second 
term of (|4Tj) one finds 

X!(( ( ^2,k 3 (7'C'2,k 4CT 'C'2,kiaC'2,k2a) ~ (Qfaa'^Mcr' ) (A^kia^teja) ) = ^ki,k 4 ^k 2 ,k 3 ~ 
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1 ^k 4 ^k 2 X k2 



^.ks+kx-ka, ^Z{C\ Ma >C 2 ^ i a'Cl M(J Ci^ 2a ) - <5 k4ikl <5 k3 ,k 2 - (48) 



Then the density-density correlation function becomes 

7V a ki.ka.ks.k* a k3 ^ kl)(i a kl)d 

( a k 3 ,d a k 4 ,d + ak 3 ,/ a k4,/)( a ki,,i a k2,rf + a ki,/ a k 2 ,/ N 



r ■> I pin l> I i" i ^k 4 ,k 3 +ki-k 2 - (49) 

[|Ok3,d| + l a k 2 ,/rJ[|ak 4 ,dr + l a k 4 ,/IJ 
In the thermodynamic limit the contributions from k 4 = k 2 and k 3 = k x in the second 
term have zero measure. Hence all k sums can be calculated without restriction. p n , n (r) is 
seen to vanish for all r whenever the k-dependence of the a kl & coefficients is negligible (for 
example in the case Itf/tfl >> 1 f° r the solution from Sec. IV. B.). It also vanishes in the 
limit |r| — > oo, where the k — > limit of the coefficients a kj 6 gives the dominant contribution 
to (jlHJ); this behavior is indeed seen in Fig. 7. 



2. Spin-spin correlation function 

Similarly, the spin-spin correlation function 

Ps,s(r) = -^E((SiS i+r ) - (S ; )(S i+r )) (50.a) 
= Ps z ,s z (r)+p s±>Sx {r), (50.b) 

is given by 



~ T/V^E E e * (kl k2)F \/- R ki-Rk2-Rk3-Rk4{( a k3,<i a k4,<i + a k 3 ,/ a k 4 ,/) X 
4iV A <J ki,--,k 4 

l5 d«k 2 ,d + a k x J a k 2 j) [((^2,k3 ( rC'2,k4 ( 7C'2,k 1 aC'2,k 2 <7) ~ (^.kaff^^) (^.kio^ka*)) 
— ((^2,k3<T C'a.kiCT ^2,ki-a^2,k 2 -a) - (^2,k3aC , 2,k4cr)(C , 2,ki- CT C'2,k 2 - ( 7)) + (51. b) 

ak 3 ,d a k 4 ,/ - a k 3 ,/ a k 4 ,<i)( a ki,/ a k 2 ,d _ a ki,d a k 2 j)(^l,k 3CT ^2,k 4CT C'2,ki CT C , l,k 2 a)}5k4,ki-k 2 +k3 



l fl k 



Ps ± ,s ± (r) = ^EK^r +r + ^r +r )-((^f>(^ +r ) + (^)(^ +r ))] (52-a) 

^A ; 

= ^E E ^ kl - k2 > v/it! kl it; k2 it! k3J R k4 {(a k3)d a k4jd + a k3i/ a w ) x (52.b) 

/iV A <T kl,-,k 4 
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( a ki,d a k 2 ,G( + a k 1 J a k 2 j)((^2,k 3CT C'2,k4-<TC'2,ki- CT C , 2,k2 CT ) - (Cl ) k. Aa C 2< \ l4 -a)(Cl j l ll - a C2^2a)) + 
(ak 3 ,dCtU,f ~ a ^i,f a U,d)( a ti,f a t 2 ,d ~ a ki,d a k 2 j)(^l,k3 CT ^2,k4-aC'2,k 1 - CT C'l,k2 CT )}^k4,k 1 -k2+k; i - 

Thus one finds ps z ,s z ( r ) — ^r,o/4 and p s± s± ( r ) = ^r,o/2 in the thermodynamic limit. 
Consequently, the spins are uncorrelated at distances |r| 7^ 0. This is a result of the 
macroscopic spin degeneracy of the ground state. 

V. ALTERNATIVE TRANSFORMATIONS OF THE HAMILTONIAN 

A. Generalized cell operators 

The transformation of the PAM Hamiltonian into positive semi-definite form in terms of 
composite operators (the unit cell operators Ai ucr in Sect. II. B. 2) is relatively independent of 
the form of A. Instead of defining linear combinations of fermionic operators inside a unit cell 
of the Bravais lattice one may also define this superposition on an arbitrary substructure of 
the underlying lattice, e.g., a cell larger than a unit cell (for example, see Fig. 8). This leads 
to more general cell operators which may then be employed to transform the Hamiltonian 
into a form similar to (J7|). The corresponding matching conditions are similar to (JSJ) but are 
now fulfilled in a different region of parameter space. In this way it is possible to construct 
and investigate various regions of parameter space. 

Below we present a different transformation of the PAM Hamiltonian based on cell op- 
erators with octahedral shape, as shown in Fig.8. Thereby it is possible to make direct 
contact with the conventional PAM 48 . This clearly shows that a transformation of the PAM 
Hamiltonian into positive semi-definite form is not linked to finite /-electron hopping am- 
plitudes in the PAM as one might have suspected from the steps performed in Sect. II. B. It 
also shows that the conventional PAM and the generalized version with extended hopping 
and hybridization amplitudes can have ground states with quite similar properties. This 
emphasizes the fact that the physical properties of an interacting electronic system need not 
depend on the precise form of the non-interacting bands. 

We consider an octahedral cell B- x at each lattice site i. The seven sites within B- Y are 
denoted by r B . = i + r aiia2)Q , 3 , where r QljQ2jQ3 = ai(l - |a 2 |)(l - \a 3 \)xx + a 2 (l - |«i|)(l - 
|«3|) x 2 + 013(1 — |«i|)(l — |a2|)x3, a T = —1,0, 1. Here x T are the primitive vectors of the 
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unit cell of the lattice. As seen from Fig.8., the seven sites r^. can be numbered by the 
indices n ai ,a 2 ,a 3 = (1 - |or 2 |)(l - [af 3 |) (aijcKij + 2) + (1 - |«i|)(l - [af 3 |) (ck 2 |ck2| + |a 2 | + 2) + 
(1 — | cti I ) (1 — I 1 ) («^3 1 1 + 4|aa| + 2) without reference to B { . The separation vector r 
between a site i and its neighbors in the cell can take 11 distinct values (x T , x r ±x T /, 2x r , 
r, t' = 1, 2, 3, r > r') over which the summation in (|2I3|) must be performed. Instead of (J1J 
the cell operator is now defined as 

7 

^Bi,a = a 6ffi,<r + a n,d^i+r Q1 , a2 , Q3 ,o-' (53) 

where n = n ai)Q , 2)Q , 3 , and the seven vectors r aii(X2)a3 = r n = r n in the sum are, in 

the increasing order of the index n, — xi, — x 2 ,xi,x 2 , — x 3 ,0,x 3 , (see Fig.8.). We note that 
a n j 7^ only on the central site (e.g., n = 6) of the cell 5i. Consequently the product 
Bg a BB i<7 does not introduce /-electron hopping terms into the Hamiltonian, implying that 
the decomposition discussed here directly applies to the conventional PAM. In addition, 
the Hamiltonian contains only on-site and nearest-neighbor hybridization amplitudes, and 
ci-electron hopping occurs between nearest and next-nearest neighbor sites. 
Instead of (J7J) the transformed Hamiltonian then becomes 

H = £ B Bi>a Bl ia + UP + E gfi , (54) 

where E g ^ = KN — UN^ — 2A r A(|a 6 j| 2 + K), K = J2n=i \o-n,d\ 2 - Furthermore, the matching 
conditions from (JHJ) transform into the following non-linear system of 19 coupled complex 
algebraic equations 

_i xx = a 6,d a 3,d + a l,d a 6,di _if: x2 = fl M a 4,rf + a 2,d a 6,d, _ ^x 3 = a l,d a 7,d + ^d^d, 
_t x 2 +xi = a 2,d a 3,d + a*l,d a i4-> _i( X2-xi - a l,d a i4 + a*2,d a l,di _t x 3 +xi = a 5,d a 3,d + «t,d a 7,d, 
_ ^x 3 -xi = a l,d a 7,d + a^M, _ ^x 3 +X2 = a *2,d a 7,d + a 5,d a 4,d, _t x 3 -x 2 = a l,d a 7,d + a l,d a 2,d, 
_ ^2xi - a l,d a 3,d> _:t 2x2 = a 2,d a 4,d 5 "^2x3 = a l,d a 7,d' -V = a^d^./j 

— = G^,d a 6,/ = — 14 2 = a2,d a 6,/ = fl6,/ a 4,d, ~ Vx 3 = a 5,d a 6,/ = a 6,/ a 7,d, (55) 

where = — V r , Ef — K — U — |a6,/| 2 , and a n> / = for n 7^ 6. The ground state 
wave function valid for N > 3N\ now has the form 

N A 
i=l 
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where, for N = 3N\, = 1 holds. Since a Ht d/a n j cannot be constant for all n — 1,2, ..7 
(see Sec. III.) \^ g ,b) describes an itinerant ground state with properties similar to those 
presented in Sec. IV. An isotropic solution is obtained for-^ 



«M = ^, a M = ^M e ^ o)) a6i/ = _ yi ^ e # (57) 

where <fi is an arbitrary phase and <pv is the phased of the hybridization amplitude Vq. 
Introducing the notation t = {tf^/tfl, v = \Vi/tf\, this solution emerges in the parameter 



U + Ef 2v 2 1 

-3t--r + ^7— T2 • (58) 



region with 



\tf\ " t ' 2tcos 2 

Since Vq and thereby Jm(Vo) is in principle arbitrary, (jSHJ) defines an entire region in the 
parameter space where the solution exists (see Fig. 9.). For Im(Vo) 7^ 0, or real Vq at 
In/^l > 6, the properties presented in Sec. IV. remain valid.— 

The results presented in this Section show that an itinerant non-Fermi liquid phase 
emerges also for real hybridization — even in the conventional PAM. 



B. Alternative decomposition using unit cell operators 

To emphasize the flexibility of our method for constructing exact ground states of the 
PAM we will now show that even if unit cell operators Aj ua are used, alternative decompo- 
sitions are possible which lead to qualitatively different ground states in other regions of the 
parameter space. For example, instead of ((7|) the Hamiltonian H defined by (jl.all.cl2j) can 
also be cast into the form 

H = E4„«A" + + K- ( 59 ) 

i,cr i 

For periodic boundary conditions the matching conditions can then be obtained from (JHJ) 
by (i) replacing J^' b ' by J^' b ' = —J b ' b> , and (ii) replacing the last two equations (for V and 
U + Ef) by Vq = ELi <, d anj and E f = K d - K f , with E'JN = K d . 

Using the new solutions of the matching conditions for arbitrary U > we will now 
identify new ground states of (jEHj) at and below quarter filling. 
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1. Insulating ground state at quarter filling 

There exist solutions of the form q n ,d/qnj = P for all n, with the additional constraint 
p = p*. In this case, the unit cell operator Aj ha contains only contributions of the form 
fi,a+pdi, a for all i and a. Introducing = [f- A - (l/p)d\ A ] + /^[/^ - (l/p)4j, where 

Hi is an arbitrary site dependent constant, one observes that 

{fj,a+pdj,aJyM} = 0, Vj,j',/iy,<X. (60) 

Since the operator Jj (fJ,i) acts only on site i, the operator P = Yii=L 0^0 does not introduce 
/-electron double occupancies. Consequently, £) i)(7 . Aj. a Ai ucr P\^) = 0, UJ2i n{+n{ J t|0) = 0, 
and the ground state becomes 

JVa I i 

i*p> = IIK/ilt - -4t) +^(/u - -4i)]i°>- ( 61 ) 

i=l P P 

The ground state (jHTj) has exactly one electron per site, and is degenerate, localized and 
globally non-magnetic. One can directly show that i?eo" TjT (0) = (see Sec. III. A. 1.), i.e., the 
state is indeed insulating. 

2. Conducting ferromagnetic ground state at quarter filling 

If q n ,d/Qn,f depends on n, A Ii)(J does not anticommute with //(//j); hence 
Y,i a A\. a A Iha P\Q)) = does not hold. In this case one may introduce a complementary 
unit cell operator— 

8 

n(a,/3,7)=l 

= (ll^la + ?2,d^i+xi,CT + ?3,d C 4+xi+X2,o- + 9 , 4,d^i+x 2 ,o- + — • + ?M^i+x 2 +x 3 ,cr) 

(62) 

with the property {Ar i)CT , = for all i, j and a, a'. The numerical coefficients are given 

by the relations 

Q*,d = wa 7j, Q*2,d = wa s,f, (lid = wa 5j, qt, d = W «6,/, 

ltd = wa 3 ,f, ql d = wa 4J} q* d = wa ljf , = wa 2 j, 

q* lf = -wa 7jd , qlj = -wa 8:d , q* AJ = -wa 5:d , q* Af = -wa 6id , 

95,/ = ~ wa 3,f, qtj = -wa 4 , d7 q? tf = -wa ltd , q* SJ = -wa 2 , d , (63) 
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where w is an arbitrary non-zero constant. Introducing Q] Ti a2 UN = Yii=i Q\. a - one nas 
J2i,(T Ai i)CT Ai ha Ql i a2 ^ N ^ |0) = for arbitrary a- x . If the sites i, i' are adjacent the oper- 
ator Q\ iCTi Qi hai introduces two electrons on the common sites of I\ and Iy. The state 
Q^ ai a2 aN |0) then provides the minimum possible eigenvalue (e.g., zero) in the presence of 
the Hubbard term only when o- x = a for all i. Consequently, the ground state becomes 

N A 

« = n<Mrl°>;*=U- (64) 

i=l 

This corresponds to a fully saturated, non-degenerate ferromagnetic phase, which is metallic 
since in double occupancies (e.g., d\+fh), empty sites, and single occupancies are 
simultaneously present.— 

Ferromagnetism in the PAM at and around quarter filling has been investigated rather 
extensively in the past already. Following variational results for Kondo lattices by Fazekas 
and Miiller-Hartmann^ ferromagnetic phases in the PAM itself were found by Dorin and 
Schlottmann^ in the limit U — > oo, and at finite U by Moller and Wolfle^ 7 - within a slave- 
boson approach. Subsequently, ferromagnetic solutions of the PAM at and near quarter 
filling were obtained within various other approximation scheme o^ 1 ^ 1 ^ 1 ^ 21 ^ 4 ^ 1 ^ . g xac j. 
results were derived in D = 1 for U = oo by Yanagisawa 67 . 

In the exact solution discussed here ferromagnetism emerges when the lower diagonalized 
band of (J2J) becomes non-dispersive (flat). We note that this can happen even if the bare 
bands of the Hamiltonian are dispersive. 

3. Ground states below quarter filling 

Exact ground states of the PAM can even be constructed for N < i.e., below quarter 
filling. For the two cases discussed in Sec.V.B.1-2. they have the same form as (jMj) and 
(jfiTjl . However, the upper limit of the products has to be replaced by N, and the additional 
geometrical degeneracy of the electrons needs to be taken into account. 

For n-dependent q n ,d/q n ,f, the ground state becomes 

W = E E E^-^kwiIK II QlJ\o) (65) 

{V N , m } {a m } m ie»iv,m 

where T>^ denotes an arbitrary domain (a subset of lattice points) of the full lattice con- 
taining iV < A^a lattice sites. A given domain V N consists of disjoint subdomains (clusters) 
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denoted by X>jv,m where m enumerates the clusters in T>n; the maximum of m is denoted by 
N m . Clearly one has T>m = Vn,i U Pjv,2 U ... U Pjv,jv TO , an d T^N,mi H T>N,m 2 — for mi 7^ m 2 . 
Furthermore, a m = ±1/2 represents a fixed, but arbitrary, spin- index in the sub-domain 
T> NtTn , and ccd n ,{t> n m },{a m } are arbitrary coefficients. From a physical point of view (J55J) 
contains disjoint, fully saturated ferromagnetic, conducting clusters of arbitrary shape and 
size, whose spin orientation is arbitrary. 

For q n ,d/c[n,f = p = p* the ground state becomes 

W = E E a^wIK II [(/^-kt)+^-^)])l°> ( 66 ) 

containing again clusters of arbitrary shape and size. But now the clusters are insulating 
and non-magnetic, containing strictly one particle per site with arbitrary spin. 

The parameter region corresponding to the ground state solutions in Sec.V.B. can be 
obtained from that derived from the matching conditions (JEJ) by the replacements t\ — > —t h r , 
V r — > —V r , U — > 0. The results are valid for all U > 0. 

VI. SUMMARY 

We presented details of an analytic scheme which allows one to construct exact ground 
states for a general class of t/iree-dimensional periodic Anderson models (PAM), including 
the conventional PAM, on regular Bravais lattices. First the Hamiltonian is cast into positive 
semi-definite form with the help of composite fermionic operators in combination with a set of 
coupled, non-linear matching conditions for the microscopic parameters of the Hamiltonian. 
Then a non-local product state of these composite operators in position space, corresponding 
to 3/4 filling is constructed which yields exact grounds in various parts of parameter space. 

Depending on the choice of the composite operators and the geometry of the building 
blocks of the lattice on which they are defined, the transformation of the Hamiltonian into 
positive semi-definite form can be performed in several ways. Thereby it is possible to 
construct exact ground states in different regions of the parameter space of the model. 

For real d, f hybridization amplitudes we constructed an insulating, non-magnetic ground 
state which is stable on several different lattice structures. Its ground state energy was shown 
to diverge at the boundary of the stability region, implying a divergence of its compressibil- 
ity. Such an anomaly is known to occur in several heavy-fermion materials. Furthermore, 
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we identified an exact metallic non-Fermi liquid ground state, characterized by one dispers- 
ing band and one upper flat band, which is stable in different regions of parameter space. 
This state is non-magnetic and has vanishing non-local spin-spin correlations in the ther- 
modynamic limit. Its density-density correlations are short ranged, and the momentum 
distributions of the electrons in the interacting ground state have no discontinuities. The 
stability regions of these ground states extend through a large region of the parameter space 
- from weak to strong on-site interactions U . 

Exact ground states with conducting and insulating properties, respectively, were also 
constructed at and below quarter filling. In particular, a conducting, fully polarized ferro- 
magnetic state was found to be the ground state at quarter filling. At lower fillings a ground 
state characterized by ferromagnetic clusters of arbitrary shape was identified. 

Our results show that ground states of the conventional PAM and of generalizations with 
extended hopping and hybridization amplitudes can have quite similar properties. 

The exact ground states discussed in this paper correspond to simple solutions of the 
coupled matching conditions for the microscopic parameters of the Hamiltonian. In view of 
their large number (e.g., 55 conditions in the case of the unit cell operators in Sec. II. B) and 
their non-linearity it almost certain that other solutions exist which then lead to yet other 
exact ground states of the three-dimensional PAM and its extensions. In view of the great 
relevance of this model for our understanding of correlated electronic systems, both on the 
level of models and real materials, more detailed investigations of the matching conditions 
derived here will be worthwhile. Finally it should be stressed that the concept behind the 
construction of exact ground states for the PAM in D = 3 presented here is quite general, 
and is also applicable to other electronic correlation models. 
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APPENDIX A: TRANSFORMATION OF THE THE PERIODIC ANDERSON 
MODEL HAMILTON! AN 



In this appendix we present details of the transformation of the PAM Hamiltonian in its 
original form with d and / operators [see (1), (2), (6)] 

i,<r r 

+ {V Q d\J- ha + H.c.) + (E f + U)n{j + UP- UN A , (Al) 
into the form with unit cell operators, (J7J), 

H = - E \A h ,° + K dN + UP- UN A , (A2) 

i,<r 

and explain how the matching conditions (JSJ) arise in this process. 
One first calculates all terms in the sum 

-Ei CT i/ iCT 4r 1)CT from (|A2|l . To identify the 
contributions generated thereby with those defined by (jAl|) one needs altogether 55 matching 
equations. This procedure will now be illustrated by four typical examples. Here we refer 
to Fig. 10 where four neighboring unit cells 1^,1^,1^, and I- J2 - each denoted according to 
the notation introduced in Sec. II.B.l - is depicted. Since the precise shape of the unit cells 
is unimportance here we use orthorhombic cells for simplicity. 

Nearest neighbor contributions: First we analyze a nearest neighbor amplitude, e.g., 
ij 2 j 5 , which appears in the kinetic energy of (jAl|) as tj 2 ,j 5 d\ 2 ^di 5 ,a- This term acts along 
the bond ( j 2 ? j 5 ) represented by the horizontal, thick dashed line in Fig. 10. In flA2|) this 
term originates from — X)i )(7 A\. a Ai ha , but only for sites i = ii , 12 5 ji , J2 in La- To calculate 
~~ Za=ii,i 2 ji j 2 ^j.^Aj. ^j one must write down the unit cell operators A\ , corresponding 
to the aforementioned four unit cells. For example one finds 

b=dj 

b=dj 

Thereby one obtains contributions of the form —c^ y d a T,dd\ 2U d- }5(T , " a t,d a 8,ddj 2fT dj 5jCT , 
-al jd a 3 , d dl 2t(T dj B ,a, and -a* id a M d} 2 ,<A^ implying fj 2h = -(a* l d a 4 , d + a* 2 , d a 3,d + al^d + 
a 6,d a 7,d)- Since (i) tf 2 j 5 is a specific example of the general hopping amplitude tf i +X2 , (ii) all 
amplitudes amplitudes have the same form (i.e., t^ 2 = tf i+ ), and (iii) this also holds for 
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the / electrons, i.e., t£ 2 = t{ i+X2 , one finds -t^ 2 = a{ b a^ b + a* 2 b a^ b + al b a S)b + a* 6b a 7)b , where 
b = dj . 

The same analysis applies to the hybridization amplitudes. Hence, for the terms 
V x2 f dl a fi+ X2 ,<r and VlfftAi+^v one finds -V*f = a* lb a^y + a* 2b a^ b , + a* 5b a 8:b/ + a* 6b a 7jb >. 

Thus we have derived the matching condition J X2 6 ' = — i X2 + (1 — ^&,&')^X2 6 '] third 
relation in (jHJ)]. 

It should be noted that the sum — J2i,a A\ a^h-p a ^ so produces the Hermitian conjugates 
of the terms represented by (jAlj) . Therefore, in addition to the result for tj 2 ■ one also finds 

*jsj2 = -( a l,d ai > d + a l,d a ^ d + a M a M + a 7,d a 6 i d), where t- 5 ^ 2 = ^ ; _ X2 = t d _^ 2 = £ X2 . 

Second-nearest neighbor contributions: We now consider a plaquette-diagonal next- 
nearest neighbor term, namely tf 2 j 5 , acting on the oblique thick dashed line in Fig. 10. 
Since the bond (12, J5) appears only in the unit cells L n and I{ 2 , only the contributions 
from — ZaAj. a Aj. ajCT in (jA2j) with i = ii,i 2 give rise to tf 2 ^ 5 d\ 2a dj 5iCT in (jAljl . The sum 
- Si=i 1 ,i a A} ucr A Iua leads to the expressions -a* 24 a 74 d\ 2 (J d^ a , and -al d a S4 d\ 2 (7 d- }5)(T , im- 
plying tf 2 j 5 = —{a* d a 8,d + a 2.d a 7,d)- Since (i) tf 2 j 5 is a specific example of the general 
hopping amplitude tfi +X3+X2 , and (ii) all the amplitudes tfi +X3+X2 have the same form (i.e., 
^n+x 3 +x 2 = ^x 3 +x 2 ) one finds — ^x 3 +x 2 = a M a 8,rf + a 2,d a 7,d- These relations also hold for the 
/ electrons. Therefore we arrive at the matching condition J X3 + X2 = a\ b a^y + a\,b a 7V [the 
sixth relation in (JHJ)]. 

Third-nearest neighbor contributions: The third neighbor contributions, located along 
the space diagonals of the unit cell, are the same in every cell. For example, studying the 
bond (ji, J5) of the cell L n) the product — A\ t a A Ih a leads to the term — a* t dCi7,ddl lja d^ 5ta . 
Hence one finds tf ljB = tf )i+X3+X2+Xl = * X3+X2+X1 = ~a*x,d a iA and, consequently, J x ' 3 + X2+Xl = 
a id a 7,d [the tenth relation in (JHJ)]. 

Single-site contributions: Since any site is common to 8 unit cells, single site con- 
tributions require the consideration of 8 neighboring unit cells. For example, the term 
VoT,i,a <?,«r/?,«T is obtained from - £ i<T A\ ua A Iua as -(En=i <,<£««,/) la,* <?L /U> implying 
=1 a n,d a nj- I n the same way we obtain the coefficient of — ^,0-^,0- as 
ELi l«n,&| 2 = #6- Therefore the terms - £ i>(T A\.^A Ih(T + # d A> in flU yield -{K d N d + 
if/iV/) + K d (N d + N f ) = (K d - Kf)N f , where N = N d + N f . Taking into account the term 
(U + Ef)N f in ((SU one obtains U + E f = K d - K f [last relation in ©]. 
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APPENDIX B: THE CURRENT OPERATOR 



In this appendix we derive the current operator as well as the charge conductivity sum 
rule for a general form of the PAM. 

The charge density operator at site i is given by p(i) = e2 ff (4Af + fi a fi,a)-> with e as 
the electron charge. The current operator j is defined as^ j = — (z/7i)[q, H], where V is the 
volume of the system, q = (1/V) J2i ip(i) is the charge polarization operator, and H is the 
Hamiltonian. With H from (jAl|) the current operator of the PAM is found as 



"' EEE E (^A--^iU r . ( B1 ) 



^ i r a b,b'=d,f 



where J^' b ' is defined in (|5|). and r is restricted to the values presented in Sec.II.B.l. We 
note that only the itinerant part of the Hamiltonian, i.e., H itm — J2 r Hiu n (r) = H — Hjj — 
ZiA( V °d\Ji,* + H.c.) + E f nQ, where 

H iUn (r) = Ei( td JlA+r,a + t f JtJi+r,a) + {V? d\Ji + r,„ + V/ d %d l+r J + H.C.} (B2) 



contributes to flBjJ) . 

Starting from the current operator in (|B1|) . the Kubo formula for the charge conductivity 

at zero temperatures becomes a T)T [uj) = iV / °° dte~ lwt ([q T (t),j T \). The sum rule for the 

charge conductivity then has the form2& J °° duRea TjT (uj) = — ^([g T) jV]) 5 where q T = q ■ 

x T /|x r | and j T — j • x T /|x T | represent the r components of q and j, respectively; here x T are 

the primitive vectors of the unit cell. For the Hamiltonian under investigation one finds 

7tp2 (r ■ x "12 ~ 
duRea TjT {u) = -7^7 E K -j-^-(H ltin {v)) . (B3) 

2kV r l X r| 

APPENDIX C: LOCALIZED SOLUTION FOR THE NON-CUBIC CASE 

In this appendix the matching conditions (|23|) from section III.B are solved for localized 
ground states in the case of non-cubic systems. We start from the observation that the 
coefficients a n ^ with n > 5 can be expressed by the unit cell diagonal hopping amplitudes 
and a n>d with n < 5 as a 5jrf = -ti 3 _ X2 _ Xl /al d , a M = -ti 3 _ X2+Xl /a^ d , a 7<d = -t X3+X2+Xl / ' a i,d, 
a 8,d = — ^x 3 +x 2 -xi/ a 2 d- With these relations and using the notations 

Afd fd. Afd +d Afd fd 

_ — ^^X 3 +X2+Xl^X 3 -X2-Xl _ X 3 +X 2 +Xj X3+X2-X1 _ X 3 +X 2 +Xj X3-X2+X1 

Tl = Td 2 ' = Td 2 ' T3 = Td 2 > 

t X2+Xl fc X 3 +X 2 ''XS+Xl 
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Afd fd Afd fd Afd fd 

„ — X3-X2-X1 X3+X2-X1 _ X3-X2-X1 X3-X2+X1 _ X 3 -X 2 +X! X3+X2 — Xj //-li \ 

T 4 = 71 2 ; T 5 = 2 ' T 6 = 13 2 'l Ui J 



leads to 



fd z iJ t d t d 

b X3-X 1 ''X3-X2 t x 2 -xi 



*X2+X1*X 3 +X 2 T?(r \ T?(„ \ T?(„\ l„ |2 _ *X2-X1*X 3 +X2 T 2 -F^) ^(je) 



= - TT ^ ^(ri)F(r 2 )F(r 4 ), |a 2 



24_ Xl r 4 ^'^y^v^ 24+xi T3F(T2) 

24+X2 ^2)^4)' ' Ml 24 3+X2 r 2 F(r 3 ) ' 



where F{x) = 1 ± yl — x. The function required for the expression of _E 9 becomes 

8 y-d 2 j-ti 2 

EL |2 _ /|_ |2 _| X3+X2+X1 \ , /i |2 | X3+X2-X1 \ 
— \\ a l,d\ H i in — J + U a 2,rf| H i in — J 



n=l 



\ a 2,d\ 

+ (|«3,d| 2 + ) + + ^^ir 1 )- ( C3 ) 

For the localized solution to be stable the quantities |a nj( i| 2 m (|C2|) . and hence the function 
F{ T i)i f|C3|) . need to be real. When -F(Tj) becomes complex the localized solution becomes 
unstable. At these points the derivative dF(x)/dx diverges. Except for accidental cancella- 
tions in (jC3|) . this implies an infinite slope of E g as a function of hopping amplitudes. 



APPENDIX D: SECTOR OF MINIMAL SPIN OF THE ITINERANT SOLUTION 

Here we present details of the calculation of spin expectation values in the itinerant case 
for minimum total spin (see Sec.IV.E.2). 
In (14 lj) the terms 

Na/2 Na 

[ n (5> iki /l)] = E /WW . /W* = E (-i)V (lll ^ + - +1 ^^ A /»), (di) 

ieD CT k {ki} P {k . } 

appear where -Fj{ k .} CT = (ft il afk i2 a---fk i / a)- Here the sum over {kj} extends over all 
sets of momentum with elements chosen from the first Brillouin zone, the sum J2p Ilr 

goes over all permutations -P{k 4 } of the momenta in each set {kj} = (ki, k^/2), and 
p represents the parity of -P{k,}- Since for i n>c G V a we have = + R, one finds 
fahl = e ^ {k * } Pv^hh where ^{ki} = R Eke{kjk. Thus, flUJ) becomes 

N A 

M = E E a^MMhk.MM) > Km,{M> = [n4 >T 4j4 {ki} /i{k i} ,jO), (D2) 
{ki} {k,} k 
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where a {ki}>{kj} = /5 {ki}iT /? {ki}a , |a {ki}i{kj} | = |a {kj}i{ki} |, and |u {ki}){k . } ) are orthogonal 
states. Using (jD2|) one finds^ 



(S x ) = (£SZ) = 0, (&) = (£§® = 0, (5 2 ) = (E^) = 0, (S) = 0. (D3) 

k k k 

To calculate expectation values of the square of the spin we use normalized wave functions 
|w{kj},{ki}) = |v{k J },{k l })/((v{k 3 },{kj|v{k 3 },{kj)) 1/2 , in terms of which the ground state can 
be written as \V g ) = £ {k . } E{k J } a '{k i },{k j }\ w {K},{*i}) > wnere a {ki},{kj} are new numerical 
coefficients. Thereby one finds 

= + S-S+)) = E{k ' } '^M^P^a - ^mm) < 

4 4E { kjE { k J }|a { k 1 },{k ;) }l 2 4 

where ^{ki}^^} is the number of common elements of the sets {kj} and {k,}. 



APPENDIX E: EXPECTATION VALUES FOR THE FLAT BAND 

In this appendix we present details of the calculation of ground state expectation values 
for the itinerant solution in Sec. IV. The ground state wave vector (|42|) is a superposition of 

states 

= n^jtniTrE^^,)]!"). ( ei ) 

k n =l V JV A kn 

where {a} = (cr^, a- l2 , a- l3 , cr^ ). By modifying the {a} sets in flElj) one obtains 2^ states, 
which obey (* 9 , W |^, {CT '}> = Det[x- hi ({X*}, {X k }, a h trj)], where Xij({U k *}, {IU k }, a i; aj) = 
5 ff .^(l/V A )E k A V^Wke^- 1 ). We see that for {</} ^ {a} the states |* 9 , M ), and |* fl>{( ,'}) 
are not necessarily orthogonal. 

1. The case X k ^ 

We first consider V k 7^ for all k. Introducing the states 

N A N h i N a i 

i*j* } > = tn^kT^jtniT^ E^ eJkni "^j]io)' (E2) 

k n =l V^va kn ^- kn 

one finds (#J ct} |# 9i{(T '}) = Det[z^{a u aj)], where Zy(cri, oj) = a: y ({l/X k }, {V k }, Oi, oj) = 
5 ffii(T ^ij, such that Det[zij{a h a\)\ = 5 {a}i{a > } . This yields (^ { g a} \^ g:{(j/} ) = S{<r},{<r'}, i-e., the 
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set {l^g^o-})} is linearly independent and for iV = 3 A/a provides a basis for H g (see (fTUjl ). 
An arbitrary state of the form of can then be written as \^ g ) = J2{a} a wl^s,w)) where 
ctjo-} are numerical coefficients^. Furthermore, one finds 

ci^Am^.w}) = e i^}i k2CT2) )' i $ Kf 2ff2) ) = tn^VAVj >< 

n=l k' 
1 1 ^Va 1 



[(-= £ ^e ik ^J )kv; )-..-(-= E ^'^ k "-^^ )( 7 =X k2 e^ i «C 2 t ikiCTi 5 (7L , CT2 ) 

V^ V A k' 11 V^ V A k / '"-i V^ V A 

1 Na I N A 

(7rE^ eik ' i,+1 4y J'-C-tttE^^^ )ll°>. ( E3 ) 

from which 

(^I^A^I^'}) = E Det[zuM,o$[ (E4) 

n=l 

follows. Here, the matrix Zij{n, of, aj) is the same as ^(of, crj), except for the matrix ele- 
ments in the n-th column which are given by z\ n = (l/NjCjiX^ /X k )e^ k21 "~ kll - ) 5 fT ;' o-,^; CT2 . 

' i ' in ' 

From (|E4J) one finds that {^/ g a "^\Cl kiU C2,]i 2 ~a\^g,{a'}) vanishes in the thermodynamic limit 
as 1/N\. In the case of o\ = 02 = cr, only {er"} = {a'} components remain in (|E4|) yielding 

X 1 Na 

(*¥ } \cl M Ate^ 9 ,w}) = Swyjrw- E <W;/ (k2_kl) S (E5) 

where X^i^o-cr; quantifies the number of a spins from {cr'}. Since J2a^^ = 1 f° r 

'in 1 in 

arbitrary cr[ n , it follows that (^J CT " } | E CT AUi A> k 2 CT l*9>'}} = S k 2 ,ki5{a"},{*>}- Using the 
notation (...) = (ty g \...\ty g ) / g \ty g ) , one finally obtains 

(E^V 1 Ak 2CT > = <w (E6) 

a 

The spin dependent expectation value {C\ ka C^M) ma y ^ e calculated, for example, by taking 
the T — > limit of (Qj kcr C2,ko-) as 

(i) = lim TriAe '^ = Urn rr ' ie '^» = SW^IAI^M), (E7) 

where A is an arbitrary operator and if 9 is defined in ()36|) . The second equality in ()E7|) 
holds since for A/a < oo excited states of H are always separated by a finite energy from the 
ground state and thus give only exponentially small corrections to the contribution of H g . 
Eq.(|K7|) therefore yields 

(QUiCTiC^kaCTa) = 2^ k 2>ki^2,Ti- (E8) 
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2. The case X k = 

Here we consider the case ak*,d = for a given vector k = k* which implies Xk* = 
(we note that the operators are then still well-defined and R^} ^ 0), with X k ^ k » ^ 0. 
Such a situation arises, for example, in the case of the itinerant solution (see Sec. IV. B) for 
large next-nearest neighbor hopping amplitude Itf/t^l < 3. Except for k* ()E6IE8|) remain 
valid the thermodynamic limit, and the ground state expectation values of the momentum 
occupation becomes 

(C'l,k<jC'l,k -) = 1, (Cl^ a C 2 ,k*a) = 0, (C f 2,k^k*<rC f 2,k^k*cr) = ^ ( E9 ) 

We note that for one vector k* the maximum total spin decreases from to — 1 

since in the sum over k n in (El, E2) the term with k* is missing. 



3. Expectation values for correlation functions 



A key problem in the calculation of correlation functions is the evaluation of expectation 
values of products of four C^ko- operators. Using (|E3|) one finds 

N A 

< ^2,k 3CT3 ( ^ 2 . k 4^4C , 2,k lCTl C'2,k 2 <7 2 |^ 9 ,{ CT '}) = S kltk J aua4 \^g,{a'}',n ^ ) + 

n=l 

N A N A 

e i^ 1 ^^^^), (eio) 

n=l m=l,mj^n 

where k 4 = k 3 + k x - k 2 . Here l^^^f 2<Ta) ) is defined in JEfl, and 



|a / [(k30- 3 ),(k 4 (T4)];[(ki(7i),(k2(T2)] 
I g,{a'},m,n 



1 N A | i iVA 

'iVA >m-l ViVA lm A/iVA 



iV A 1 iVA 



A k / 



N a 



■m + l 



A k / 



A k / 



>n+l 



TrE^^^ )]|o> 



(Ell) 



contains only contribution with n ^ m. The expectation value of the density-density cor- 
relation function can be calculated from (|E10Jl . with o\ = o 2 = <J a , 03 = 04 = o&, where a 
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summation J2cr a Ha b has to be included. Introducing the notations Z^'^ = Z^S^^, where 
-^k3,ki = ^k4^k 2 /(^k 3 ^ki^A), and using the procedure leading to (jE6|) one then finds 

( Y ^2,k 3CTf) C , 2,k4(7, ) C , 2,ki ( 7 a C'2,k2<7a) = <5ki,k 4 £k2,k 3 + #ki,kA 3 ,k 4 - ^k 3 ,ki^k 4 ,k 3 +ki-k2- (E12) 

The same expression is obtained if we calculate the expectation value in the T — > limit, 
as described in connection with (jE7l - IE8|) . 

In the case of the spin-spin correlation function for the S z components one again has 
°\ — a 2 — cr a and o" 3 = o" 4 = o"6, but the sum over the spin indices must be separately 
performed for crj, = a a = a and = — a a = a. Using (jE7HE8IE"T0j) . one finds in the T — > 
limit 

1 1 -k k 

(Y ^2^C2,k 4 aC 2 ,ki^ 2 ^) = ^ki,k 4 #k 2 ,k 3 + 2^ki,k 2 ^k 3 ,k 4 - 2^k3 4 ,ki^k4,k 3 +ki-k 2) 

(52 ALsa^^aAUi-aCVka-a) = ^ki^kg.k* - ^^ka.ki^.kg+ki-ka- (E13) 

To calculate the spin-spin correlation functions for the S x , S y components, (|E10J) must 
be evaluated for cr 3 = cr 2 = —&i = — cr 4 = cr, and the summation must be performed. In 
the T — > limit one finds 

(5Z ( ^2,k 3( TC'2,k4-aC'2,k 1 -aC'2,k 2 <T) = ^kj ^ka.kj, , 
a 

(X!^lk 3 aC'2,k4- CT C'2,k 1 -aC'l,k2<T)^k4,k 3 +k 1 -k 2 = ^ka.ks^ki.kj- (E14) 
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FIG. 1: A unit cell I\ connected to an arbitrary site i showing the primitive vectors x T and indices 
n of the sites in /. Arrows depict some of the hopping and hybridization amplitudes (J = t, V) 
defined within I\. 
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FIG. 2: Ground state energy per lattice site. (|25|) . in units of U as a function of y = (t^Ail f° r 
the localized solution in a simple cubic crystal for |i{/if| = 0.01. As seen from the plot, E g /N^ is 
finite at y c = 1/2 but has infinite slope. 
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FIG. 3: Surface in parameter space representing the stability region of the conducting ground state 
discussed in Sec. IV. A. 




1 2 

FIG. 4: Tilted unit cell located at site i discussed in Sec.IV.B. Numbers represent the intra-cell 
numbering of lattice sites. Only those bonds are presented along which hopping of d electrons 
occurs. 
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FIG. 5: Surface in parameter space representing the stability region of the itinerant solution derived 
in the case of the conventional PAM, For — * Oj the surface asymptotically approaches 

the (\VL/4\,(U + Ef)/\4\) plane. 
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FIG. 6: Momentum distribution functions n^ CT , n^ a for the conducting solution discussed in 
Section IV.B, for = 5.0, = 0.5, tf > 0, and k 2 = h = 0, where k T = kx r . The plot 

presents the behavior in the first Brillouin zone for k\ £ [— vr,7r]. For other k directions a similar 
behavior is found. 
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FIG. 7: Density-density correlation function for |tf/#j| = 4, |Vi/if| = 0.25, tf > for the 
itinerant case obtained from (|49j) in the thermodynamic limit. The nine-dimensional integration 
was performed by a Monte Carlo method using 69 points. The distance r was taken in r = 1 
direction, and is expressed in units of the lattice constant a. 
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n=4 



n=5 



FIG. 8: Octahedral cell defined at lattice site i as discussed in Section V.A. The r = 1,2,3 axes are 
represented by arrows with broken lines, x T are indicated by thick full line arrows, and n represents 
the cell independent notation of sites inside the octahedron. 
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FIG. 9: Surfaces in parameter space above which the conducting phase discussed in Section V.A. 
is stable, a) |Vi|/|if| < 0.5; (U + Ef)/\tf\ is seen to diverge as l^l/l^il approaches zero, b) 
l^i l/l^i I > 0.5. In contrast to Fig. 5 where V\ is imaginary, V\ is real here. 
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FIG. 10: Four neighboring unit cells in D = 3. Each cell (e.g., I\) is denoted by one of the lattice 
sites at which it is located (e.g.,i); see Sec. II.B.l. For simplicity, orthorhombic unit cells are used. 
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